SLAC-PUB-13650 
SU-ITP-09/22 



New Fixed- Target Experiments to Search for Dark Gauge Forces 

James D. Bjorken,^ Rouven Essig,^ Philip Schuster,^ and Natalia Toro^ 

^Theory Group, SLAG National Accelerator Laboratory, Menlo Park, GA 94025 
^Theory Group, Stanford University, Stanford, GA 94305 
(Dated: June 3, 2009) 

Fixed-target experiments are ideally suited for discovering new MeV-GeV mass (7(1) gauge bosons 
through their kinetic mixing with the photon. In this paper, we identify the production and decay 
properties of new light gauge bosons that dictate fixed-target search strategies. We summarize 
existing limits and suggest five new experimental approaches that we anticipate can cover most of the 
natural parameter space, using currently operating GeV-energy beams and well-established detection 
methods. Such experiments are particularly timely in light of recent terrestrial and astrophysical 
anomalies (PAMELA, FERMI, DAMA/LIBRA, etc.) consistent with dark matter charged under a 
new gauge force. 



I. NEW GAUGE FORCES 

The interactions of ordinary matter establish that 
three gauge forces survive to low energies. Two strik- 
ing features of these forces — electroweak symmetry- 
breaking at a scale far below the Planck scale and ap- 
parent unification assuming low-energy supersymmetry 
— have driven model-building for a quarter-century. 
But the strong and electroweak forces need not be the 
only ones propagating at long distances. Additional 
forces, under which ordinary matter is neutral, would 
have gone largely unnoticed because gauge symmetry 
prohibits renormalizable interactions between Standard 
Model fermions and the other "dark" gauge bosons or 
matter charged under them. 

There is an important exception to the above claim: 
new "dark" Abelian forces can couple to Standard 
Model hypercharge through the kinetic mixing operator 
IF^F'^'^, where F'^^ = d[^A'^-^ and A' is the dark gauge 
field H]. If the A' is massive. Standard Model matter ac- 
quires milli-charges proportional to e under the massive 
A'. Kinetic mixing with e ~ 10~^— 10~^ can be generated 
at any scale by loops of heavy fields charged under both 
U{iy and U{1)y, and the A' can acquire mass through a 
technicolor or Higgs mechanism. A mass scale near but 
beneath the weak scale is particularly well-motived — 
U{iy symmetry-breaking may be protected by the same 
physics that stabilizes the electroweak hierarchy |i2j. In- 
deed, if the largest symmetry-breaking effects arise from 
weak-scale supersymmetry breaking, then the 11(1)' sym- 
metry breaking scale is naturally suppressed by a loop 
factor or by y/e, leading to MeV to GeV-scale A' masses 

HIEISISJIB]. 

An A' can be produced in collisions of charged particles 
with nuclei and can decay to electrons or muons. The 
production cross-section (cr^') and decay length (7cr), 

CTA' 100pb(e/10~'')^(100MeV/mA')^ (1) 
7CT ~ 1 mm(7/10) (lO'^/e)^(100MeV/mA') (2) 



vary by ten orders of magnitude for the e's and masses 
niA' we consider. This wide range calls for multiple ex- 
perimental approaches, with different strategies for con- 
fronting backgrounds. Beam-dump searches from the 
1980's exclude the low-mass and small-e parameter range, 
and other data constrains large e. In this paper we sug- 
gest five scenarios for fixed-target experiments sensitive 
to distinct but overlapping regions of parameter space 
(see Figure [T]). Together they can probe six decades in 
A' coupling and three decades in A' mass with existing 
beam energies and intensities. 

Dark matter interpretations of recent astrophysical 
and terrestrial anomalies provide a further impetus to 
search for new t/(l)'s. Annihilation of dark matter 
charged under a new 11(1)' into the A' can explain the 
electron and/or positron excesses observed by PAMELA 
[7], ATIC [8], FERMI [9], and HESS [11 HI] (see 
e.g. [11 m m 1151 iMl HZl HSl Hg). If the dark matter 
is also charged under a non-abelian group, then its spec- 
trum naturally implements an inelastic dark matter sce- 
nario 20J , thereby explaining the annual modulation sig- 
nal reported by DAMA/LIBRA [21 , 22 and reconciling it 
with the null results of other experiments [H [501 HSl HI] • 

In view of these suggestive data and the abundant the- 
oretical speculation surrounding them, insight from new 
experiments is clearly called for. New probes of weakly 
mixed MeV-GeV U{l)'s directly probe the low-energy 
structure of these scenarios, where the nature of their in- 
teractions is most manifest. As such, the experiments we 
advocate here are complementary to upcoming gamma- 
ray observations (see e.g. [5S]|25]) and to the next gener- 
ation of direct detection experiments [571 1251 1^ [Ml [5T| 
that will shed light on the scattering of dark matter. 



Direct Tests of Low-Mass Gauge Sectors 

Constraints on new A's and the reach of different ex- 
periments are summarized in Figure [l] To begin, low- 
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FIG. 1: Left: Existing constraints on an A' . Shown are constraints from electron and muon anomalous magnetic moment 
measurements, Oe and a^, the BaBar search for T(35') "ffi'^fi^, three beam dump experiments, E137, E141, and E774, 
and supernova cooling (SN). These constraints are discussed further in Section III Right: Existing constraints are shown in 
gray, while the various lines — light green (upper) solid, red short-dashed, purple dotted, blue long-dashed, and dark green 
(lower) solid — show estimates of the regions that can be explored with the experimental scenarios discussed in Section IV A - 
|IVE[ respectively. The discussion in |IV| foc uses on the five points labeled "A" through "E" . The orange stripe denotes the 
"D-term" region introduced in section jil A[ in which simple models of dark matter interacting with the A' can explain the 
annual modulation signal reported by DAMA/LIBRA. Along the thin black line, the A' proper lifetime cr — SO/im, which is 
approximately the r proper lifetime. 



energy e^e~ colliders are a powerful laboratory for the 
study of an A' with e > 10~^ and mass above ^ 200 
MeV, particularly in sectors with multiple light states 
[321 IMl iMl [351 136^- Their reach in e is hmited by lu- 
minosity and irreducible backgrounds. However, an A' 
can also be produced through bremsstrahlung off an elec- 
tron beam incident on a fixed target '345 . This approach 
has several virtues over coUiding-beam searches: much 
larger luminosities, of 0(1 ab~^/day) can be achieved, 
scattering cross-sections are enhanced by nuclear charge 
coherence, and the resulting boosted final states can be 
observed with compact special-purpose detectors. 

Past electron "beam-dump" experiments, in which a 
detector looks for decay products of rare penetrating par- 
ticles behind a stopped electron beam, constrain > 10 
cm vertex displacements and e > lO"*". The thick shield 
needed to stop beam products limits these experiments to 
long decay lengths, so thinner targets are needed to probe 
shorter displacements (larger e and m^' ) . However, beam 
products easily escape thin targets and constitute a chal- 
lenging background in downstream detectors. 

The five benchmark points labeled "A" through "E" 
in Figure [l] (right) require different approaches to these 
challenges, discussed in Section [IV| We have estimated 
the reach of each scenario, summarized in Figure [l] 
(right), in the context of electron beams with 1-6 GeV 
energies, nA-/xA average beam currents, and run times 
~ 10^ s. Such beams can be found for example at the 



Thomas Jefferson National Accelerator Facility (JLab), 
the SLAG National Accelerator Laboratory, the electron 
accelerator ELSA, and the Mainzer Mikrotron (MAMI). 

The scenarios for points A and E use 100 MeV-1 GeV 
electron beam dumps, with more complete event recon- 
struction or higher-current beams than previous dump 
experiments. Low-mass, high-e regions (e.g. B and G) 
produce boosted A' and forward decay products with 
mm-cm displaced vertices. Our approaches exploit very 
forward silicon-strip tracking to identify these vertices, 
while maintaining reasonable occupancy — a limiting 
factor. At still higher e, no displaced vertices are re- 
solvable and one must take full advantage of the kine- 
matic properties of the signal and background processes, 
including the recoiling electron, using either the forward 
geometries of B and G or a wider-angle spectrometer (e.g. 
for point D). Spectrometers operating at various labora- 
tories appear capable of probing this final region. 

We focus on the case where the A' decays directly to 
Standard Model fermions, but the past experiments and 
proposed scenarios are also sensitive (with different ex- 
clusions) if the A' decays to lighter [/(I) '-charged scalars, 
and to direct production of axion-like states. 

Outline 

In Section [ll| we summarize the properties of A' pro- 
duction through bremsstrahlung in fixed-target coUi- 
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sions. Constraints from past experiments and from neu- 



trino emission by SN 1987A are presented in Section III 



In Section [TVl we describe the five new experimental sce- 
narios and estimate the Hmiting backgrounds. We con- 
clude in Section [V] with a summary of the prospects for 
new experiments. More detailed formulas, which we use 
to calculate our expected search reaches, and a more de- 
tailed discussion of some of the backgrounds, are given 
in Appendices |A) 15] and [O . 



II. 



THE PHYSICS OF NEW U(l) VECTORS IN 
FIXED TARGET COLLISIONS 



Theoretical Preliminaries 



Consider the Lagrangian 



C 



\F''^'^Fl^ + m:'^,J^^^^, (3) 



where £sm is the Standard Model Lagrangian, = 
d^^A'^Y and A' is the gauge field of a massive dark U{1)' 
gauge group [T]. The second term in ([3| is the kinetic 
mixing operator, and e ~ 10~® — 10^^ is naturally gen- 
erated by loops at any mass scale of heavy fields charged 
under both C/(l)' and J7(l)y; the lower end of this range 
is obtained if one or both C/(l)'s are contained in grand- 
unified (GUT) groups, since then e is only generated by 
two-or three- loop GUT-breaking effects. 

A simple way of analyzing the low-energy effects of the 
A' is to treat kinetic mixing as an insertion of p^gi^^—p^Pi, 
in Feynman diagrams, making it clear that the A' couples 
to the electromagnetic current of the Standard Model 
through the photon. This picture also clarifies, for ex- 
ample, that new interactions induced by kinetic mixing 
must involve a massive A' propagator, and that effects 
of mixing with the Z-boson are further suppressed by 
Equivalently, one can redefine the photon field 
A^ A^+eA'f^ as in ^7\, which removes the kinetic mix- 
ing term and generates a coupling eA^ J^j^ D eeA'^J^^^ 
of the new gauge boson to electrically charged particles 
(here e = eycos^v^)- Note that this does not induce 
electromagnetic millicharges for particles charged under 
the A'. The parameters of concern in this paper are e 
and niA'- 

We now explain the orange stripe in Figure [l] — see 
[31 m [S] for more details. In a supersymmetric theory, 
the kinetic mixing operator induces a mixing between 
the D-terms associated with J7(l)' and U{1)y- The hy- 
percharge D-term gets a vacuum expectation value from 
electroweak symmetry breaking and induces a weak-scale 
effective Fayet-Iliopoulos term for 11(1)'. Consequently, 
the Standard Model vacuum can break the C/(l)' in the 
presence of light [/(I) '-charged degrees of freedom, giving 
the A' a mass, 



niA' 



(4) 




FIG. 2: A' production by bremsstrahlung off an incoming 
electron scattering off protons in a target with atomic number 
Z. 
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FIG. 3: (a) 7* and (b) Bethe-Heitler trident reactions that 
comprise the primary QED background to A' £'^£~ search 
channels. 



where gn, gY, and 32 are the the U{1)' , U{1)y, and 
Standard Model SU{2)l gauge couplings, respectively, 
and mw is the W-boson mass. Equation Q relates 
e and m^' as indicated by the orange stripe in Figure 
[1] for gD ~ 0.1 — 1. This region is not only theoret- 
ically appealing, but also roughly corresponds to the 
region in which the annual modulation signal observed 
by DAMA/LIBRA can be explained by dark matter, 
charged under the ?7(1)', scattering inelastically off nuclei 
through A' exchange. We therefore include these lines for 
reference in our plots. 



B. A' Production in Fixed- Target Collisions 

A' particles are generated in electron collisions on a 
fixed target by a process analogous to ordinary pho- 
ton bremsstrahlung, see Figure [2] This can be reli- 
ably estimated in the Weizsacker- Williams approxima- 
tion (see Appendix [A] for more details) |38l [39l l40] . 
When the incoming electron has energy Eq, the differ- 
ential cross-section to produce an A' of mass niA' with 
energy Ea' = xEq is 



da 

dxd cos 9 A' 



IP 



(1-2;+-)- 



a;2 x(\~x)m\,{Elxe\:) 



C/2 



(5) 



where Z is the atomic number of the target atoms, 
a ~ 1/137, Qa' is the angle in the lab frame between the 
emitted Al and the incoming electron, the Cog (~ 5 — 10 
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for niA' 500 MeV) depends on kinematics, atomic 
screening, and nuclear size effects (see Appendix |A] and 
Figure 10 therein), and 



Uix,9A') = E^x9%+m\, 



+ mx 



(6) 



is the virtuahty of the intermediate electron in initial- 
state bremsstrahlung. The above results are valid for 



1/ < 1. 



(7) 



Dropping me and performing the angular integral, we 
find 



da 
dx 



3(1-2:) 



Cog 



(8) 



of massless 



5l with m\, = 



The x-scaling and singularity structure 

bremsstrahlung |3D] is recovered from 
(the polynomial factor differs because of finite m^' cor- 
rections to the matrix element) , but differs from the mas- 
sive A'-strahlung in several important ways. We empha- 
size that these properties are not particular to any matrix 
element, but rather is a kinematic property common to 
all heavy-particle emission: 

Rate: For most x, U{x,0) ~ iti\,, so that the total A' 

3 2 

production rate is controlled by Therefore, it 

A' 

is suppressed relative to photon bremsstrahlung by 

C 9 
A' 

Energy: A' bremsstrahlung is sharply peaked at a; w 
1, where [/(x, 0) is minimized. When an A' is 
produced, it carries nearly the entire beam en- 
ergy — in fact the median value of (1 — x) is 



Angle: A' emission is dominated at angles 0a' such that 
U{x,9a') ^ 2U{x,0) (beyond this point, wide- 
angle emission 

falls as 1/6*^, )• For X near its median 
value, the cutoff emission angle is 



max 



y/rnA' 



m„ m 



3/2' 



En 



E, 



3/2 



(9) 



which is parametrically smaller than the opening 
angle of the A' decay products, ~ itia'/Eq. The 
approximation of coUinear emission is justified in 
many calculations. 

Because these features apply to bremsstrahlung of any 
massive boson, there is a simple approximate equiva- 
lence between constraints on the A' and constraints on 
an axion with mass rua — ruA' and decay constant 
/o ~ 1.7 X 10~^GeVCe/e, at least when the constraints 
come from coupling to electrons. Here the coupling con- 
stant of axions to electrons is assumed to be Carrie/ fai 
where Ce is a model-dependent coefficient. 



The total number of A' produced when Ng electrons of 
initial energy Eq scatter in a target of T radiation lengths 
is 



dN_ NoXo 
dx ~ ^ A 



da 



dtI{Ei;Eo,t) 



X E, 



dx' 



c'=E^,/Ei 



(10) 



where x' = Ea' / Ei , Xq is the radiation length of the 
target, A^o — 6 x lO^'^mole"^ is Avogadro's number, A 
is the target atomic mass in g/mole, and / is the energy 
distribution of electrons after passing through t radiation 
lengths. 

After the A' is produced in the target, it will travel 
for some distance before decaying back into Standard 
Model particles (we will assume throughout this paper 
that no other decay channels into particles charged un- 
der the U{1)' are available). The proper lifetime of the 
A' is 



CT 



1 _ 3 

r NcsiriA'ae'^ 
1-4 ^ 2 



80 iim 



10- 



100 MeV \ 
rriA' ) 



(11) 



^, where R is defined 



where we have neglected phase-space corrections and AoS 
counts the number of available decay products (A'off = 1 
for VTLA' ^ 2m^ when only A' e+e~ decays are possi- 
ble, and 2 -I- R{mA') for uia' > 2m 

to be the energy dependent ratio °'^^^^.|, 
A' decays will thus create displaced vertices behind the 
target. While cr determines the typical impact parame- 
ter for these displaced tracks, their vertex displacements 
are controlled by (for the typical kinematics with x ^ 1), 



hadrons) \ 



io 



'SEi 



0.8cm / Eo \ 
N,s UOGeV/ 



10- 



^^100MeV\^ 



ruA' 



(12) 



where we have again neglected phase-space corrections. 



C. Approximate Total Rate Formulas 



From equations ([s]) and (12 1, we can obtain simple 
approximate expressions for the rate of A' production 
in scattering off thin targets (with T <C 1) and thick 
"dump" targets (T 3> 1). These crude approximations 
are only correct within about one order of magnitude, 
but they are useful in quickly mapping out regions in 
the large logarithmic parameter space. In our results we 
use more accurate expressions that also include detector 
acceptances, as presented in Appendices \K\ and [B| 

In the thin-target limit T <C 1, the beam is not sig- 
nificantly degraded as it passes through the target, and 
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I{Ei,Eo,t) « S{Ei — Eq). In this case the total A' pro- 
duction rate scales as 



A 



T 



Cog = N^CT€^ 



(13) 



where C ~ 5 is only logarithmically dependent on the 
choice of nucleus (at least in the range of masses where 
the form-factor is only slowly varying) and on m^/ , be- 
cause, roughly, Xq oc (see Appendix [b| and [41j). For 
example, for a Coulomb of incident electrons 



N 
C 



10'' 



T 



0.1/ VlO 



2 /lOOMeVy 
V mA' J 



(14) 



For a thick target (T ^ 1), production is dominated near 
the front of the target and 



N NX' < 



2 



(15) 



with C w 10. When the typical lifetime Iq exceeds the 
length L to the detector, a fraction ~ i/^o decay before 
the detector, and the number of A! observed is indepen- 
dent of TYlA''- 



No 



bs 



NX'e 



NX'c 



El 



(16) 



Note that multiple interactions in the target degrade the 
beam energy significantly and induce A' transverse mo- 
menta ~ 10 mrad (GeV/i?o)^ for A' production in the 
first radiation length. These transverse momenta can be 
significant for low-energy dumps. 

For subsequent discussions, it is useful to translate the 
signal yields into rates as a function of beam and target 
parameters. A process X with cross section (y{X) occurs 
with a rate 



$(X) - 0.7 MHz 



'E ' Ibeam 

nA 



a(X) 



(17) 



where Ibeam is the average current, and T ^ 1 is the tar- 
get thickness in units of radiation lengths. For example, 
ruA' = 100 MeV and e = lO"* gives aA' O.OIZ^ pb, or 



a rate of « 0.007 Hz 



In contrast, Bethe- 



Heitler pair production (Figure |3| has a total cross 
section ai+i- ~ Z^fih, or rate « 0.7 MHz '^'^j^'" 

for a ~ 1 GeV electron beam. Bethe-Heitler pair produc- 
tion is thus a potential background for any A' search, and 
so our experimental scenarios will be strongly infiuenced 
by the need to remove them. 

In appendix [C] we discuss the kinematics of Bethe- 
Heitler production relative to A' production in some de- 
tail, and sketch out a set of selection cuts that can be 
used to suppress the otherwise prohibitively large Bethe- 
Heitler backgrounds (Figure [3](b)). However, a minimal 
contribution to the background is obtained by replacing 
the A' by a 7* (Figure [sja)). To see this, we re-insert 



the dilepton invariant mass into the fully differential 
cross-section and consider integrating over a mass win- 
dow Sm, with r <C Sm <^ m. The A' and 7* matrix 
elements are related by the substitution 



1 



— (mj^, + iniA'Ty 



(18) 



All other terms in the cross-section are slowly varying in 
this window. Treating them as constant, the integration 
over the mass window bounded by m ± ^ is straight- 



forward. Substituting (111, we obtain the ratio of fully 



differential cross sections for A' to 7* production in this 
mass window, which is also an upper bound on the total 
signal to background. 



dcr{X A'Y l+l-Y) 
da{X -> 7*r ^ l+l-Y) 



2Nja 



v om 



(19) 



where N f is the number of available decay species for 
the A' , and 5ra is the width assigned to the 7* process. 
Equation (19 1 summarizes the maximum achievable sig- 



nal to background ratio that any experiment can achieve 
in an A' i'^i^ search using only kinematics, with the 
decay vertex unresolved. 



III. BEAM DUMP CONSTRAINTS AND 
SENSITIVITY OF CURRENT EXPERIMENTS 

In this section, we discuss existing constraints on the 
e versus niA' parameter space, which are summarized in 
Figure [T] 

For TTiA' > 2m^, a search for T(3S') — > jA' 7^+/i^ 
by the BaBar collaboration rules out e > 10^'^ (see 
also [5?), while the electron and muon anomalous mag- 
netic moments rule out the low-mass-high-e region [43j . 

Strong constraints are also obtained from electron 
beam-dump experiments searching for MeV-mass axions. 
The strongest constraints come from the El 37 [H] and 
E141 [45] experiments at SLAC, and the E774 [46] ex- 
periment at Fermilab: 

SLAC E137 dumped 30 C of electrons at 20 GeV into 
aluminum targets [Hj. Beam products traveled 
through a 200 m hill and an additional 200 m of 
open region before hitting a (3 m)^ detector. No 
candidate events were observed. The contour in 
Figure [1] represents an expected signal of 10 events 
(we have idealized the detector as a circle of radius 
1.5 m). 

SLAC E141 dumped 2 x 10^^ electrons at 9 GeV into 
a 12-cm tungsten target, with a 10-cm tungsten 
target used for calibration [15]. The detector was 
located 35 m from the dump, and the analysis re- 
quired observing a single decay product carrying 
over 0.5 times the beam energy with angular ac- 
ceptance set by a 7.5-cm pipe. Based on the back- 
ground rates reported by the experiment, the ex- 
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elusion in Figure [T] represents an expected signal of 
1000 events. 

Fermilab E774 dumped 0.52 x 10^° electrons at 275 
GeV onto two 28-radiation-length-thick (about 
19.6-cm) stacks of tungsten plates 05]. The overall 
target length, including veto counters behind the 
target, was 30-cni. An electromagnetic calorime- 
ter with an angular acceptance of about 20-cm was 
placed 7.25 m downstream from the dump. The 
trigger required an energy deposition of at least 
27.5 GeV and no signal from the veto counters. 
Based on the results reported by the experiment. 
Figure [T] represents an expected signal of 17 events. 



The approximate formulas given in Section |II C| are suf- 
ficient to understand the shape and magnitude of the 
beam-dump limits shown in Figure [ij they are bounded 
above by a diagonal along which many A' may be pro- 
duced, but all decay within the shielding that stops the 
beam, and from below by a line of diminishing rate, which 
is diagonal if the typical decay occurs before the detec- 
tor position, and approximately horizontal if the average 
decay length £o is larger than the length scale L of the 
experiment (see equation 16 1). Similar limits can be de- 
rived for alternative A' decay modes, for example if the 
A' decays to dark-sector higgses with typical proper life- 
time that scales as rather than e~^. 

Supernova cooling places a significant constraint on 
lower e and lighter niA' ■ A proper accounting of super- 
nova limits on the A' is beyond the scope of this paper, 
but we outline a simple estimate based on scaling simi- 
lar results for axions jlT]. The hot core of a collapsing 
supernova can cool through production of A' if they de- 
cay > 10 km from the point of production (the mean 
free path is typically longer than the lifetime). However, 
neutrino observations of SN1987A confirmed an energy 
loss over 5-10 seconds of 1 — 4 x 10^^ erg. Following [47], 
we require the energy loss in A' emission not to exceed 
Iq53 erg/s. We take the A' luminosity per unit energy 
from the core to be 



dL 

dEA 



1 

TsN 



(6 X 10'" erg/s)e''e^ 



(20) 



for Ea < TsN = 30 MeV, which is suppressed by 
TsN/fnp (where nip is the proton mass) relative to the 
axion rate |471 I48j because the vector emission matrix 
element is proportional to w^, whereas the axion emis- 
sion matrix element approaches a constant as u — > 0. 
We impose an additional Boltzmann suppression e^^-^^"^ 
for Ea > Tsn and multiply by the fraction /{Ea) = 
g-io km/^o ^]^gj^ leave the supernova core. Requiring that 
the total luminosity not exceed 10^'^ erg/s, we can ex- 
clude the lower-most region in Figure [T] We emphasize 
that the luminosity obtained by scaling is only correct 
within an order of magnitude. An error in the cross- 
section would affect the lower limit in e proportionally, 
but the upper limit only logarithmically. 



Constraints from other experiments are all contained 
within the limits from the experiments discussed above. 
For example, the region constrained by the SLAG search 
for milli-charged particles, which also used an electron 
beam, is contained within E137 [49, . Proton beam dumps 
can produce A' in radiation directly from the proton or in 
radiation from electrons produced by the nuclear shower. 
Both processes produce A' of much lower energy than 
the primary proton (hard bremsstrahlung off the proton 
is suppressed by the proton's finite size, and the shower 
electrons are quite soft). These softer A' typically decay 
inside the dump. Therefore, proton dumps such as the 
GHARM experiment at GERN ^ do not exclude new 
regions, though they do overlap significantly with the 
E137 exclusion. Likewise, experiments dumping proton 
beams for other purposes (e.g. neutrino experiments such 
as MINOS and MINIBOONE) have Uttle or no potential 
reach beyond E137. 

We have also considered potential limits from A' 
production off cosmic rays impinging on Super-K, 
AMANDA, and IGE-GUBE detectors, which is domi- 
nated by bremsstrahlung off muons near ground-level, 
which must only survive 1 km to reach the detectors. 
The potential sensitivities of these experiments are con- 
tained within the E137 excluded region. 



IV. SCENARIOS FOR NEW EXPERIMENTS 

The parameter space that new experiments must cover 
spans a huge range. The A' production cross section and 
decay width vary as and thus vary over ten orders 
of magnitude. Our purpose in this section is to explore 
experimental scenarios appropriate to different parame- 
ter ranges. For the sake of definiteness, we organize the 
discussion around parameter points labeled "A" through 
"E" in Figure [l] Each choice suggests a different experi- 
mental approach, described in the appropriately labeled 
subsections. 

We do not intend here to provide detailed designs; this 
is the task of those who would actually do the experi- 
ments. We do attempt to show that the exploration of 
this parameter space is experimentally feasible and offer 
some guidance regarding how to choose design param- 
eters in order to optimize the experimental sensitivity. 
Because the electron beams at Jefferson Laboratory ap- 
pear to be an attractive choice for such experiments, we 
have been guided in our considerations by the beam spec- 
ifications available there. However, we expect that other 
attractive options exist elsewhere. 

There is a natural dividing line in the parameter space 
(cf. Figure [iJ, corresponding to an A' proper lifetime 
CT ~ 80 /im, comparable to that of the r lepton. Longer 
lifetimes allow in principle the determination of a sepa- 
rated decay vertex, while much shorter lifetimes do not. 
Since determination of the detached vertex is a strong 
experimental signature, the experimental techniques nat- 
urally differ in the two regimes. Beam-dump searches. 
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including our first (and last) scenario, are appropriate to 
much longer lifetimes. The region near the dividing line 
has not yet been explored, and micro-vertex detectors 
appear quite promising in this range. The second and 
third scenarios we describe assume this technique. 

For very short lifetimes of the A' , the experimental sig- 
nature is identical to electromagnetic trident production 
e + Z^3eoTe + Z^e + 2iJ, (where Z is the target), 
shown in Figure [Sj The simple upper bound discussed 
in equation ( [T9| in Section |ll on the ratio of the A' fully 
differential cross-section to the background trident rate 
implies that high statistics and resolution are required to 
have any chance of observing the A' . Fortunately, such 
a regime does appear to overlap with the capabilities of 
Jefferson Laboratory spectrometers. Moreover, the up- 
per bound is attained in a sizeable region of the differ- 
ential phase space where A' production is dominant, if 
appropriate kinematic cuts are applied on the final state 
leptons (see Appendix |C| . Therefore this parameter re- 
gion seems in principle to be accessible. Our fourth ex- 
perimental scenario deals with this case. We will also 
show that the second and third scenarios we discuss have 
some interesting new reach in parameter space when used 
at lower luminosity as high-resolution forward spectrom- 
eters. 

For the smallest values of e, the primary considera- 
tion is simply producing enough A's to study experimen- 
tally, so beam-dump experiments are the technique of 
choice. However, it becomes very challenging to design 
beam dumps with average power exceeding 1 megawatt 
(MW). There appears to be a small window of oppor- 
tunity available for such a search, which would increase 
the reach beyond that of E137. This comprises the fifth 
scenario that we discuss. 



Low Power, 10 cm Tungsten Beam Dump; 

e = 10"'^; ruA' = 50 MeV 



We consider a 200 MeV primary electron beam inci- 
dent on a 10 cm tungsten target. Downstream (beyond 
a thin, but dense shielding wall, if necessary), is an in- 
strumented decay volume containing a combination of 
tracking planes, electromagnetic calorimetry, and scintil- 
lator triggers. With the chosen values of e and toa', the 
laboratory decay length of a typical A' of momentum 160 
MeV is about 5 cm. The produced A's are contained in 
an angular cone of order 125 mrad. Therefore the track- 
ing system in an experimental region no more than 40 
cm (see Figure [4]) downstream of the front of the dump 
need have transverse dimensions no more than 10 cm to 
identify the A' decay vertices and measure the angles. 
Since the decay angles of the electron/positron pair are 
of order 250 mrad, the calorimeter transverse dimensions 
can be very modest. 

In this scenario, the total yield per incident electron of 
A's containing at least 80 percent of the beam momen- 
tum is about 9 x 10"^^ per electron dumped. If the front 



of the fiducial decay volume can be located immediately 
behind the 10 cm target, 5% of these A's decay outside 
the target for a rate of 4 x 10^^^ per electron dumped. 
If a thicker shield is necessary to stop soft photons, the 
yield remains large: 0.1% of the A' decay beyond 30 cm 
from the front of the dump, for one observable A' decay 
per 5 X 10^^ electrons dumped. In the conservative con- 
figuration, with a total of 30 cm of material, a yield of 
order 30 events would be observed for 0.3 coulombs of 
electrons dumped (300 nA in an experiment of duration 
10^ seconds). This requires a modest 60 watts of beam 
power on the tungsten dump. 

The length of the fiducial decay region need only be 
20 centimeters to capture the majority of the A's emerg- 
ing from the dump into the decay region. The compact 
nature of this decay volume suggests the possible use of 
silicon strip detectors for the tracking system. 

The question of backgrounds must of course be ad- 
dressed. A fast, dense tracking system seems to be ap- 
propriate. With a readout rate at least 10 MHz, and 
with a continuous incident beam such as exists at Jef- 
ferson Laboratory, there would be about 30000 electrons 
dumped per readout cycle. The shower products should 
be absorbed efficiently, and it should be very rare that a 
prompt calorimeter signal of more than 100 MeV energy 
deposition occurs, especially because the time resolution 
of a scintillator trigger/electromagnetic calorimeter sys- 
tem will be much better than 100 ns. The residual prob- 
lems, beyond the scope of this sketch, probably have to 
do with neutrons and soft photons or x-rays. 

Evidently, if we lower niA' and increase e in such a 
way as to decrease neither the rest-frame decay length 
nor the production rate, the experiment will be easier, 
since jct is larger and more A' will decay in the detector 
volume. In Figure |4] we present our rough estimate of the 
region of parameter space accessible to the experiment as 
described. 

To extend the reach to larger masses and smaller e, 
higher beam energies are required. Once the threshold 
for electro-production of muons and hadrons has been 
crossed, the experiment may require a thicker shield, and 
much higher energies appear advantageous. For example, 
consider raising ttia' and lowering e by a factor of 2, to 
100 MeV and 5 x 10~^. For a beam energy of 6 GeV, 
the decay length is 2 m. If a 3 m decay volume can 
be positioned with its upstream end within 4 m of the 
front of the dump, then a larger fraction of the produced 
A's can be detected than our original example. This in 
turn lessens the beam-intensity requirement; an average 
current of 100 nA, which leads to a dump power of under 
1 kilowatt, appears to suffice. 

The detector geometry can simply be a longitudinally 
stretched version of the previous case, with transverse 
dimensions again quite small, of order 15 - 20 cm. How- 
ever, new backgrounds appear. Muons will penetrate the 
decay volume as well as electromagnetic showers initi- 
ated within the hadronic cascade. One leading candidate 
for background troubles comes from electro-production 
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FIG. 4: Left: Experimental scenario for benchmark point A (e ~ 10~^, m^/ ~ 50 MeV). An electron beam is incident on a 10 
cm thick tungsten target. Behind the target is a 10 cm (or thicker) shield followed by an instrumented decay region consisting 
of a combination of tracking planes, electromagnetic calorimetry and scintillator triggers. Right: Reaches of the high- and 
low-energy dump configurations described in Section |IV A[ delineated by regions with 10 or more events and the following 
configurations — Blue (inner) Solid Contour: 0.3 C total charge dumped with a 200 MeV electron beam, a 20 cm shield, and 
a detector with 5 cm radius 50 cm behind the front of the target. The lepton pair must have total energy exceeding 100 MeV. 
Blue (inner) Dashed Contour: same configuration, but with no shield. Creen (outer) Solid Contour: 0.1 C (100 nA beam x 
10^ s) total charge dumped with a 6 GeV electron beam, a 3.9 m shield, and a detector with 10 cm radius 7 m downstream. 
The lepton pair must have total energy exceeding 3 GeV. Creen (outer) Dashed Contour: same configuration, but with 0.9 m 
of shielding. Cray contours and Orange Stripe: exclusions from past experiments (E137 and E141) and the region that explains 
DAMA/LIBRA in a simple model — see Figure [l]for more details. 



of the p, with a leading charged pion from the rho de- 
cay undergoing a charge-exchange reaction into a tt'' a 
few radiation lengths in front of the detector region. We 
have used the experience obtained in E141 to make rough 
estimates, which indicate that such backgrounds are sur- 
mountable. But the soft backgrounds such as neutrons 
and hard x-rays also need to be carefully studied. 

B. Thin Target and Double Arm Spectrometer; 

e = 3 X 10"^; m^' = 200 MeV 

Modern micro-vertex detectors allow much better life- 
time resolution than the above example. When e is in- 
creased from the previous example, the rate of A' pro- 
duction per incident electron increases, and a thin target 
can be used instead of a beam dump. For the parame- 
ters of interest here, we consider a 0.1 radiation length 
tungsten target. We choose a 6 GeV beam with an av- 
erage current of 100 nA. Downstream of the target is 
a two-arm mini-spectrometer with silicon strip detectors 
as the tracking elements, backed up with fast calorime- 
ter/scintillator triggers. 

With these parameters, the A' production rate (before 
acceptance) out of the target is about 10 per hour. The 
angular divergence of the A' beam is only about 5 mrad. 
The laboratory decay length is about 1 cm, and the de- 
cay products of the A' have an average angle of about 



35 mrad from the beam axis. A spectrometer with polar 
angle coverage of 20 to 55 mrad and 50% azimuthal an- 
gle coverage has about 25% acceptance for the A' decay 
products. The trigger requirement includes the demand 
that the energies in each of the calorimeters are between 
1 and 5 GeV, with the sum between 5 and 6 GeV. The 
tracking system must identify one track in each arm that 
points to the calorimeter hit (if the calorimeter is seg- 
mented) and is consistent with a decay-vertex origin. Af- 
ter reconstruction, additional kinematic constraints pro- 
vide rejection power. In Figure [5] we show the reach 
of this experimental scenario for various geometries and 
different beam currents. 

A major background is simultaneous elastic coulomb 
scattering in each arm. An elastically scattered electron 
deposits 6 GeV in the calorimeter, and is rejected, but 
the singles rate must be below one per timing window 
(100 MHz or less for fast calorimeters). This require- 
ment is safely met by the beam intensity quoted above. 
The elastic-scattering radiative tails will contribute to 
the trigger, but at a significantly lower rate of 10 kHz or 
so. Other sources for background triggers, such as Bethe- 
Heitler pair production (cf . Figure |3| , lead to smaller or 
comparable trigger rates. When one of the two scattered 
electrons scatters again in the first layer of silicon, the 
intersection of the two reconstructed tracks is displaced. 
We find that the rate for these fake vertices is adequately 
suppressed if the first layer is placed close to the target. 
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FIG. 5: Left: Experimental scenario for a small two-arm spectrometer for benchmark point B (e ~ 3 x 10^^, m^/ ~ 200 
MeV). An electron beam is incident upon a thin 0.1 radiation length tungsten target. A small two-arm spectrometer with 
silicon-strip trackers and a fast calorimeter or scintillator trigger is downstream from the target. Signal events are identified by 
requiring a displaced vertex ~ 1 cm behind the target. More details are given in the text. Right: Regions corresponding to 10 
or more events within acceptance in 10^ sec for three different geometries. From right to left: 6 GeV electron beam at 100 nA 
(0.1 C delivered), with angular acceptance from 20 to 55 mrad and aim long detector (solid red line); 6 GeV beam at 5 nA 
(5 X 10^^ C delivered), with angular acceptance from 10 to 27 mrad in a 2 m-long detector region (dashed darker red line); and 
2 GeV beam at 0.5 nA (5 x 10"'' C delivered) with the same geometry as the dashed red line (solid dark red line). In all cases, 
we require that the A' carry at least 83% of the beam energy, the track impact parameters at the target exceed 50 fim, and 
the reconstructed vertex displacement exceed 1 cm. We assume 50% (j> coverage. Gray contours and Orange Stripe: exclusions 
from past experiments (E137 and E141) and the region that explains DAMA/LIBRA in a simple model — see Figure [l]for 
more details. 



within 5 — 10 cm. 

Another basic requirement is that the occupancy in the 
tracking system be acceptably low. High-resolution sili- 
con strip detectors are beneficial in this regard. Within 
a cone of opening angle of 10 mrad at a distance of 50 
cm downstream of the target, we estimate that the den- 
sity of electrons and photons produced in the target with 
energy above 1 MeV is of order 10^/cm^/s |58]. In this 
scenario, the silicon is placed further from the beam, but 
this rate serves as a rough upper bound, which would give 
one percent occupancy for a 1 cm x 25 /im strip. While 
these numbers are encouraging, a serious simulation is 
certainly required. 



C. Silicon Strip Layers in a Diffuse Electron Beam; 

e = lO""*; ruA' = 50 MeV 

At even higher e and lower masses, there exists the 
option of halving the number of silicon strip tracking ele- 
ments and placing them directly into a defocused primary 
electron beam of low intensity. For this study, we choose 
the beam size to be about 1 cm x 1 cm and the beam 
energy to be 1 GeV. The beam intensity is hmited by 
silicon occupancy to about 10^ e~/s, if we require occu- 
pancy of about 1% in 1 cm x 25 /im strips with a timing 
window of 20 - 50 ns. 



Triggering is again accomplished by a calorimeter, with 
a strategy similar to case B and the same limitations. For 
A' masses of 20-50 MeV, decay opening angles ~ 20 — 50 
mrad are anticipated, so the calorimeter must extend 
close to the beam. For simplicity we consider an an- 
nular calorimeter with angular coverage above 20 mrad 
(for example, located at 2.5 meters from the target, with 
inner radius of 5 cm). The beam electrons emerge from 
a 0.1 radiation- length tungsten target in a Moliere dis- 
tribution, with typical transverse momenta of 5 MeV. 
Therefore less than 1% of the electron beam hits the 
calorimeter, leading to a 1 MHz singles rate, which 
is high but manageable for a trigger requiring two hits. 

With these parameters the A' production rate is about 
1 every ten hours. Off-line track reconstruction can 
be used to remove the backgrounds associated with 
the Coulomb scattering pile-up and other background 
sources, in particular Bethe-Heitler pair production from 
the target. The quality of the experiment will depend 
crucially on the precision of the vertex reconstruction 
using the silicon strip information. Our sample point 
has typical impact parameter ~ 160 /Ltm and laboratory 
decay lengths of order 2.3 mm, which should be cleanly 
resolvable. The sensitivity of this configuration, assum- 
ing several different resolutions, is illustrated in Figure 

El 

For smaller masses, the calorimeter must be placed at 
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FIG. 6: Left: Experimental scenario for benchmark point C (e ~ 10 m^' ~ 50 MeV). Silicon strip tracking elements, 
together with a 0.1 radiation length (300/im) tungsten target directly behind one of the elements, are inserted into a 1 GeV 
diffuse (1 cm x 1 cm) electron beam of intensity ^ 10* e~/s. Triggering is accomplished by an annular calorimeter with 
angular coverage above 20 mrad (e.g. 2 cm inner radius, 1 m downstream) by demanding three coincident hits carrying the 
beam energy. Signal events give rise to measurable impact parameters for the leading two tracks, and the excellent tracking 
provided by this design exploits this feature to reject background. Invariant mass reconstruction can provide an additional 
search variable (see Sec. IV D I. More details are given in the text. Right: Concentric purple contours: Regions with detectable 
signal yield > 10 events, background rejection of ~ 10~^ (yielding S/ B ^1), and an impact parameter of at least 33 /im, 66/im, 
or 150/im, respectively, for the contours from the outside in. We assume a run time of 10® s at 10* e~/s. Red Dotted Contour: 
Analogous sensitivity with lower average current (10^ e~/s) and a smaller calorimeter aperture (10 mrad). Thin black dashed 
line: a rough estimate of the total region of sensitivity that could be accessible to this geometry using both displaced-vertex 
discrimination and invariant mass search windows with good momentum resolution (see Sec IV D[ ). Gray contours and Orange 
Stripe: exclusions from past experiments (E137, E141, E774, electron and muon anomalous magnetic moments, and T(3S') 
resonance searches) and the region that explains DAM A/LIBRA in a simple model — see Figure [l]for more details. 



a narrower angle or the beam energy reduced. In cither 
case, the Moliere scattering becomes more acute. On the 
tails of the Moliere distribution, one can compensate by 
lowering the intensity of the beam. At low beam inten- 
sities, a fast scintillator/calorimeter trigger system will 
resolve the passage of individual electrons in the beam 
(in a CW machine like CEBAF). Therefore, if the scintil- 
lator/calorimeter system is segmented (e.g. scintillating 
fiber calorimetry) , the trigger requirement can be simul- 
taneous deposition of the beam energy in more than one 
detection element — typically three. For larger masses, 
the beam intensity would have to be increased, and the 
silicon-strip occupancy presents a sharp barrier. 



D. High Resolution, High Rate Trident 



Spectrometer: e = 3 x 10 



= 1 GeV 



Large A' masses present two challenges: a low produc- 
tion rate and short A' lifetime. In the absence of a dis- 
placed vertex, the A' can only be observed as a small peak 
on the electromagnetic trident background. Reducing 
these backgrounds as much as possible is essential here. 
Additionally, targets with somewhat lower Z than tung- 



sten are preferable in this high A' mass range in order 
to maintain charge coherence in scattering. For definite- 
ness, we shall discuss the di-muon final state, though it is 
arguable that the electron-positron final state is prefer- 
able. 

As discussed in Section |TT] the trident background 
arises from two subprocesses, which we call radiative and 
Bethe-Heitler (c.f. Figure |3|. The radiative process gives 
an upper bound on the ratio of signal to background as 
in equation (19 1. The Bethe-Heitler process has a much 
larger (~ lOOx) cross-section than the radiative trident 
process due to coUinear logarithmic enhancements in the 
6^67 splitting and sub-process 77 — > /i/i. These en- 
hancements can be avoided by demanding kinematically 
symmetric ^ji decay products carrying the majority of 
the beam energy, and by demanding that the recoiling 
electron (if it can be identified) scatter at a wide angle. 
This preserves the large logarithm in the forward-peaked 
A' production cross-section, while regulating all logs in 
the Bethe-Heitler process. These selections are discussed 
further in Appendix [C] 

In addition to the trident processes, radiation of real 
photons by incident electrons, and their subsequent con- 
version in the target must be considered. This process 
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FIG. 7: Left: Schematic diagram of an experimental scenario for benchmark point D (e ~ 3 x 10~*, m^' ~ 1 GeV). An electron 
beam with an energy of ~ 6 GeV and a current of about 100 /iA - 200 /xA is incident upon a 0.1 radiation length aluminum 
target. A wide-angle high-resolution spectrometer allows triggering on events in which one electron and one positron carry 
most of the beam energy. The signal is distinguished from background events with the help of various kinematic selection cuts 
(relatively symmetric final state and possible recoil electron tagging) and a "bump hunt" — see text and appendix [c] for 
further details. Right: Various estimates of the possible reaches of a wide-angle spectrometer, with (bottom) and without 
(top) tagging vertices displaced by > 1 cm to reject background. In each case, the outer thin black line represents a significant 
total rate, with no geometric acceptance requirements (S/\/B > 5 in the no-vertex (top) region, 10 or more events in the vertex 
(bottom) region). The thick blue curve shows the reach when decays are required to land more than 200 mrad away from the 
beam line, and the inner dotted curves assume an additional 1% signal efficiency from acceptance. In these two cases, each 
curve represents the total reach obtained by running at several beam energies. Gray contours and Orange Stnpe: exclusions 
from past experiments (E137, E141, E774, electron and muon anomalous magnetic moments, and T(35') resonance searches) 
and the region that explains DAM A/LIBRA in a simple model — see Figure [llfor more details. 



is naively enhanced by 0{T/a) relative to Bethe-Heitler 
trident production, but can be rejected effectively with 
the same kinematic cuts. It is, of course, reducible by 
thinning the target, which allows a compensating in- 
crease in average beam current. We have not consid- 
ered pile-up processes, but assume they are small when 
the three products are required to reproduce the beam 
energy within resolution. 

For this scenario, we consider a 0.1-radiation- length 
aluminum target in a 4 GeV beam. The total yield of A's 
is roughly 10"^^ per incident electron. If we assume an 
average beam current of 250 /j,A (beam power of 1 MW) 
and an experimental duration of 10^ sec, the total rate of 
A' production is of order one per second, or > 10^ per 
experiment. These are emitted in a cone of size ~ 100 
mrad, with decay products at opening angles near 250 
mrad and the recoiling electron at a rather wide angle, 
0.5 radians. The yield of background tridents having 
a di-muon mass within one percent of the A' mass is. 



according to (19 1, about 300 times larger, or 3 x 10^ per 



experiment. The estimated cumulative sensitivity of this 
configuration, and similar ones obtained by lowering the 
beam energy down to ~ IGeV, is illustrated in Figure 
|7] To obtain the contours in this figure, we require that 
S/VB > 5, i.e. {S/ebBo) x S > 25, where S is the 
signal rate, and B is the background rate, Bq, times the 



background rejection efficiency eb- We use equation (19 1 
to obtain S/Bq, and choose reasonable values for e^. 

The signal rate above is, indeed, larger than necessary 
for the A' resonance to be statistically significant. A less 
ambitious (and perhaps more realistic) experiment would 
also suffice for discovery. There are at least three ways 
to back off from this scenario. One way is evidently to 
improve the mass resolution. A second way is to reduce 
the beam intensity, keeping the acceptance complete. A 
reduction in beam current by a factor of 100 would still 
leave a viable signal. The third way is to reduce the ac- 
ceptance; a one percent acceptance by itself would again 
leave a viable signal. 

Optimization involves a choice of a combination of 
these factors. Jefferson Laboratory looks like an espe- 
cially appropriate venue for this scenario, with two spec- 
trometers with very good electron momentum resolution. 
In particular, the small-acceptance, high-rate spectrom- 
eters in Hall A has momentum resolution of order 10^"* 
and the large-acceptance Hall B CLAS detector has elec- 
tron momentum resolution better than 1% |51j . There- 
fore it would seem that using an electron-positron pair for 
the A' decay products may make more sense than using 
a di-muon pair. However, we feel further investigation is 
best done with the aid of expertise within the Jefferson 
Laboratory experimental community. 
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Again, it will be important to detect the wider- angle re- 
coil electron in order to reduce backgrounds. This re- 
quires increasing the tracking and calorimeter coverage 
to around 250 mrad for the two-arm spectrometer B as 
well as for the diffuse-beam, coUinear configuration C. 

The gap between the spectrometer and resolved vertex 
regimes coincides with A' decay lengths of order the la 
vertex resolution of a given detector. In this regime, one 
can gain sensitivity by imposing a loose impact param- 
eter/vertex requirement that reduces background rates 
by a few orders of magnitude. Clearly, this strategy is 
plagued by the difficulties of both the spectrometry and 
vertexing approaches, and would only be attempted in a 
later stage when the detector and backgrounds are well 
understood. We do not discuss it further. 



mA'/GeV 



FIG. 8: Estimated reach for the geometries of scenarios B 
(thick dashed red) and C (thick solid and long-dashed pur- 
ple) in "bump-hunt" operation with upgraded mass resolu- 
tion. For reference, the approximate expected reach of clean 
vertex-based searches discussed in sections B and C is given 
by the thin, black dashed contour and the thin black solid 
contour, respectively. Thick Dashed Red Contour: The com- 
bined reach (with S/ \fB > 5) of the two two-arm spectrome- 
ter geometries introduced in Section [IV B| and Figure [S] with 
a 6 GeV beam. Thick Solid Purple Contour: S/vB > 5 
for the combined reach of the different "diffuse beam" sce- 
narios in Section |IV C| assuming 1% mass resolution. Long 
Dashed, Thinner Purple Contour: Same as the thick solid 
purple contour, but using vertexing to reduce backgrounds, 
with an assumed "optimistic" rejection of 10~^ for an impact 
parameter cut of 10 /xm applied to the lepton pair. Gray con- 
tours and Orange Stripe: exclusions from past experiments 
(E137, E141, E774, electron and muon anomalous magnetic 
moments, and T{3S) resonance searches) and the region that 
explains DAMA/LIBRA in a simple model — see Figure [l] 
for more details. 



High-Resolution Forward Spectrometers at Lower 
Luminosity 

As we lower mA' while keeping e large, the production 
angles may become too low for big spectrometers such as 
those found in Hall A at JLab. The geometries of the de- 



tectors sketched in Sections IV B and IV C may be more 
appropriate. The growth of the production cross-section 
as 1/ m^, compensates for the lower average currents de- 
manded by these scenarios (10-100 nA for the two-arm 



spectrometer B (described in Section IV B I and 10 e /s 



when tracking planes are inserted in the beam as in C 
(described in Section IV Cl). 

In scenarios B and C, in the interest of simplicity, we 
did not assume a high mass resolution for the spectrom- 
eters. However, a high mass resolution, e.g. via addition 
of magnetic fields, is not at all impractical in principle. 
Assuming a 1% mass resolution, the potential sensitiv- 



E. High power, Low Energy Beam Dump; 

e = 5 X 10~®; tua' = 50 MeV 



Values of e below the E137 limit require very intense 
beams simply to produce enough A's to detect. Beam 
power limitations force one downward in beam energy. 
We choose a 200 MeV beam of electrons with an average 
current of 5 mA, representing a beam power of 1 MW. 
The A' production rate is about 2 x 10^^^ per electron 
dumped. The laboratory decay length is about 2.5 km. 
The divergence of the A' beam is about 100 mrad. We 
consider a decay region 5 m long, with its front end lo- 
cated 5 meters downstream of the dump. A tracking 
system, perhaps in the style of E137, with transverse 
dimensions 2 m x 2 m is distributed throughout the de- 
cay volume to capture the decay products from the A' 
decays. It is surrounded by electromagnetic calorimetry 
designed to efficiently capture the electrons and positrons 
emergent from the tracking volume (see Figure [9]) . With 
these parameters, the yield of detected A's is marginal 
— about 5-10 per 10^ seconds. 

This scenario has not been optimized, and other ver- 
sions less awkward can be contemplated. However, the 
region of the exclusion plot that is covered by any such 
experiment will be modest. Therefore a real design is 
likely to be opportunistic. If the experiment can be run 
parasitically for a long period of time, the benefit-to-cost 
ratio may rise sufficiently high to make such an effort 
attractive. 

If one contemplates utilizing a higher energy dump of 
1 MW power, then the number of electrons dumped de- 
creases and, keeping uia' and e fixed, the yield of A's 
decreases as well. The region of sensitivity at A' masses 
of about 200 MeV should, to a fair approximation, merge 
with the reach of E137, as estimated in Figure [9j 
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FIG. 9: Left: Schematic diagram of a beam dump design for benchmark point E (e ~ 5 x 10~*, m^/ ~ 50 MeV). A 200 MeV 
electron beam with a large current of about 5 mA (delivering 1 megawatt in power) is incident upon a thick tungsten target 
that together with shielding is about 5 m in length. Behind the shielding is a decay region 5 m long, consisting of a tracking 
system (2 m x 2 m transverse to the beam line) and surrounded by electromagnetic calorimeters — see text for further details. 
Right: Solid Red Contour: 10 events with A' energies above 100 MeV after the experiment has run for 10^ s (5000 C total 
charge dumped). Gray contours and Orange Stripe: exclusions from past experiments (E137 and SN1987A) and the region 
that explains DAMA/LIBRA in a simple model — see Figure IT] for more details. 



V. DISCUSSION AND CONCLUSIONS 

In this paper, we have described five scenarios for 
fixed-target experiments that probe kinetically mixed 
C/(l)'s with MeV-GeV masses. Kinetic mixing of size 
e ~ 10~^ — 10~^ between a fight gauge boson and the 
photon can be generated by loops of particles at any mass 
scale, with the magnitude determined by the structure of 
high-scale physics. An MeV-GeV mass for the A' can in 
turn be generated from the weak scale, especially in a su- 
persymmetric context. An A' in this mass range is also of 
interest as a possible explanation of several current dark 
matter anomalies. 

The parameter space of A' mass m^/ and mixing e 
has been constrained from two corners by existing data. 
Beam dump experiments and supernovas exclude the 
low-mass, small-e region. Larger e ~ 10""^ — 10~^ are con- 
strained for a broad range of masses by lepton anomalous 
magnetic moments and B-factory searches. 

The five approaches we have described cover the re- 
maining parameter space using fixed-target experiments 
of various geometries and 200 MeV-6 GeV beams. A 
natural extension of past beam dumps, with modest in- 
tensity and 10 cm-1 m length, can fill in the crevice of pa- 
rameter space between past beam dumps. Beam dumps 
are not well suited to searching for A' with less displaced 
decays. For these parameter ranges, thin-target experi- 
ments are required. 

Any thin-target experiment must contend with the 



backgrounds from electromagnetic electron scattering 
and trident production, which can be tackled with a 
combination of kinematics and displaced vertex selec- 
tion. Depending on m^/, more forward or wide-angle 
geometries are called for, and small-scale silicon micro- 
strip tracking can be utilized to isolate displaced de- 
cays. We have considered three such scenarios: a forward 
two-arm spectrometer, a coUinear detector in a diffuse, 
low-intensity beam, and a wide-angle spectrometer. To- 
gether, they are sensitive in the range e ^ 10^^ — 10^^ 
for A' masses from 10 MeV to 1 GeV. 

The wide-angle scenario is of particular interest, be- 
cause existing spectrometers can cover a large fraction of 
its reach. The Hall A spectrometers and the CLAS de- 
tector [ST] at JLab seem well suited for initial searches, 
and other labs may have comparable capabilities. 

Searches at low e, below the reach of the dump ex- 
periment E137, are limited by practical rate limitations. 
Power above a megawatt (MW) is difficult to sustain, 
making e ^ 10~^ — 10~^ inaccessible with beams of any 
energy in under a year of running. Our fifth scenario 
saturates this limit, with a 200 MeV MW dump, which 
can possibly be accommodated at the JLab Free-Electron 
Laser accelerator. 

When combined with existing limits, these five scenar- 
ios can either confirm the existence of new U{1) gauge 
forces at low masses or close the door on their most likely 
parameter range. 

For masses below the electron threshold, very different 
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experimental techniques are called for. These have been 
developed in [52l [53] . 

We have restricted our discussion to the simplest sce- 
nario: a single U{1) gauge boson that decays directly to 
electrons. In a larger "dark sector" , somewhere between 
this minimal scenario and the full complexity of Standard 
Model physics, decays within the dark sector dominate 
[31] . These dark-sector cascades can return some or all of 
the A' energy to Standard Model-charged particles, with 
lifetimes controlled by for vector bosons and much 
longer lifetimes cx e^* for scalars. The limits and reaches 
discussed here apply directly to any spin-1 bosons in the 
dark sector that decay directly to a lepton pair. These 
experiments are also sensitive to dark sector cascades in- 
volving spin-0 states, with appropriately deformed ex- 
clusion regions not discussed here. Besides frameworks 
with kinetically mixed U{1), these experiments are sen- 
sitive to direct production of light (pseudo) scalars (e.g. 
[Ml [5S1 [Sni [57]). It is likely that related designs more 
optimally cover these scenarios. 

We have focused here on experimental approaches tai- 
lored to A' searches in electron beams, but analyses in 
this spirit may be possible with existing data, for ex- 
ample by using beam-halo impacts in collider experi- 
ments or neutrino production beams and detectors such 
as those at Fermilab and KEK. We also have not ex- 
plored the potential of muon-beam experiments, which 
may be ideal for searches for ^'s above the muon mass, 
which are produced with rates comparable to those for 
an electron beam, but with much lower electromagnetic 
backgrounds. 



significantly alters both the kinematics and the rate of 
the process. The qualitative behavior has already been 
summarized in Section HI! 

We want to calculate the ^'-production cross-section 



dajeip) + ZjP,) ^ ejp') + A'jk) + ZjPf)) 
dE A' d cos 9 A' 



(Al) 



in the Weizsacker- Williams approximation following [33 
[Ml 130] , where k — {Ea' , k) is the momentum of the out- 
going A' , 9 A' is the angle of its momentum relative to 
the incoming electron momentum p in the lab frame, 
p — {Eq,p) and Pi — (Mi,0) are the initial momenta 
of the electron and the target of mass Mi and atomic 
number Z, and p' — (E'^jf) and Pf are the outgoing 
four- momenta of the electron and target, which are inte- 
grated over. 

In the frame of the incoming electron, the rapidly mov- 
ing atom sources a cloud of effective photons, off which 
the electron scatters to radiate an A'. Though these 
photons are spacelike, their virtuality is small compared 
to other invariants in the problem (for example mA'), 
so that the interaction of the electron with the target 
is dominated by transverse polarizations. Therefore, it 
is related to the cross-section for real-photon scattering, 
e{ph{q) e{p')A'{k) with q = Pi - Pf by [40 



dajp + P, ^ p' + k + Pf) ^ /ax\ f Epx^A' 
dEA'dcos9A' \ TT J \{l — x) 

da{p + q ^ p' + k) 



d{p ■ k) 



(A2) 



where 
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APPENDIX A: A' PRODUCTION FORMULAS 

In this appendix, we first present the cross-section for 
the production of the massive U{1)' "dark photon". A' , 
by initial- or final-state radiation off a single electron 
hitting a fixed target of atomic number Z. This process 
is analogous to photon bremsstrahlung, except that the 
coupling of the A' to electrons is e • e and the A' mass, 
mA', is much larger than the electron mass, mg, which 



X = Ea'/Eq, 
t ^ ~a\ 



(A3) 



We specify the kinematics at t = tmin and the effective 
photon flux ^ below. Note that t is not one of the 
Mandelstam variables for the 2^2 process, which will 
be denoted by — see below. 

For a given A! four-momentum fc, the virtuality t has 
its minimum value tmin when gis coUinear with the three- 
vector k — p. Solving the mass-shell conditions p'^ = 
{q + p- k)^ = ml and Pf = {Pi - q)"^ = Mf with the 
coUinear geometry, and keeping only leading effects in 



El 



9a' 



M 

E'' 



(with \q\ defined below), we find 



U 



''mm — Imin 



U 



2Eo{l - x) 

2 



2Eoil-x)J ' 



(A4) 



, (A5) 



(A6) 



where 



U = U(x,9a') = E^9l,x + mi,- — - + mix. (A7) 

X 
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At this kinematics, 



— u — — U2 = 2p ■ k — m\, = U, 



S2 — ml — 2p' ■ k + m\, = 



U 



l-x 



t2 = ip-p') 



f\2 



Ux 

l-X 



(A8) 
, (A9) 

(AlO) 



where S2, ^2, and U2 are the Mandelstam variables for the 
2^2 process. The cross-section for the 2^2 process 
is therefore 



du 2 '^^ 
d{p ■ k) dt2 



57r(s2 — m%)-^ 



32 



—u 
s 



2m^,<2 



\ + {\^xf + 



2(1 - xYm] 
IP 



Ux ^ 
\ — X / 



(All) 



where we have dropped the ^-dependence of ^ and terms approximation to the cross-section |A1| is given by 
of order in |A^p. Therefore, the Weizsacker-Wilhams 



da?. 



EqX dx dcos 9 A' 



f/2 



(1 



C/4 



Ux 



(A12) 



where /3a' = -y/l — m?^, /El. The ^-differential cross 



section is obtained by integrating (A12l with respect to 
9 A' (we will see below that x actually depends on 9a', 
but this can be neglected to excellent approximation). 
The first term in square brackets integrates to 

(A13) 



U{x, 



In the limit ttia' — * 0, this becomes the standard pho- 
ton bremsstrahlung result with a ^ singularity, while the 
second term in the square brackets vanishes. However, 
finite ttia' regulates this singularity, and in the case of 
interest, namely ttia' ^ me, we have U{x, 0) 



.2 1- 



The second term integrates to 
that 



6U{x,0) 



+ 0{mlU^), so 



d(J^-^2 
dx 



(Sa^e^;^/?^,) 
X (^1 - x - 



H- m^x 



(A14) 



This has an approximate soft electron singularity, regu- 

2 

lated by the electron mass at (1 — x)ci — Though 

A' 

not explicit in this formula, our approximations also 
break down if the electron energy (\ — x)Efj < \q\] this 
also regulates the cross-section, cutting off log(l — x) at 

(1 — x)c2 = "w- Since one cutoff or the other is always 



larger than their geometric mean m^/EQ, the A' is always 
produced from a relativistic electron. The x-integrated 
cross-section is therefore 



i ah^(3A' 



3 



log ^ 




(l-x 

















(A15) 
(A16) 



As we have noted in Section |IT] the characteristic angle 
of A' emission is set by U{x,9a') — C/(x, 0) ~ C/(x,0), 



so 9 A' ~ — , where the median value of 1 — x is 

1 — X ~ max (y^il^, T^) ■ This is parametrically smaller 

than the angle of the A' decay products with respect to 
the incoming electron, namely ~ mA'/Eo. 

We turn next to the definition of x, which is an effective 
flux of photons integrated from t — tmin to tmax, the 
total center-of-mass energy of the collision. We refer the 
reader to [351 for more details. 

For a general electric form factor G2(t), 



dt- 



-G2{t). 



(A17) 



(We note that the other form factor, Gi{t), contributes 
only a negligible amount in all cases of interest.) Al- 
though the virtual photon propagator squared, 1/i^, is 
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dominated a,t t = tmin, the final-state phase space is 
proportional to dt(t — tmin), so that virtual photons at 
all scales contribute to A' production. As discussed in 
[38 , 39 , the physical upper bound may be set not by the 
centcr-of-mass energy, but by tmax ~ , at which the 
full 2^3 matrix element begins to shut off. 

For most energies in question, G2{t) is dominated by 
an elastic component 



aH 



1 



l + t/d 



(A18) 



where the first term parametrizes electron screening (the 
elastic atomic form factor) with a — 111 Z^^^^ /me, and 
the second finite nuclear size (the elastic nuclear form 
factor) with d = 0.164 GeV^^-^/s. We have multiplied 
together the simple parametrizations used for each in 
|38j . The logarithm from integrating (A17) is large for 
t-min < d, which is true for most of the range of inter- 
est. However, for heavy A' , the elastic contribution is 
suppressed and is comparable to an inelastic term. 



G2,in{'t) 



a'H 



1 + a'2i 



(1 



Z, (A19) 



0.71 GcV^ . 



where the first term parametrizes the inelastic atomic 
form factor and the second the inelastic nuclear form 
factor, and where a' = 7Ti Z~'^/^ /mf,, rup is the pro- 
ton mass, and Hp = 2.79 |38j . This expression is valid 
when t/inip is small, which is the case for ttia' in the 
range of interest in this paper. One can show that the 
contribution from the other inelastic nuclear form factor 
Gi{t) is neghgible. 

At high masses, these simple parameterizations of the 
form factors are uncertain at the order-of-magnitude 
level. Using G2,ei + G2,m in (|A17|), and setting tmin = 
(TO^,/2£'n)^ tmax = ml 
ureHni 



we obtain xl ^ shown in Fig- 



APPENDIX B: PRODUCTION FROM THIN AND 
THICK TARGETS OF A S WITH FINITE 
LIFETIME 

Consider an electron beam with energy E'q incident on 
a target. The total number of A' that get produced in the 
target with energy Ea' = xEq and decay at a distance z 
behind the front edge of the target is given by 



dN 
dxdz 



Ne^^ I dEi I dtI{Ei;Eo,t) 

IE., JO 



A 

E{) da 



dP{z ~ ^ i) 
d^^ 



where Eq is the incident energy, the number of inci- 
dent electrons, A'o = 6.02 x 10^'^ mole~^, p and Xq are the 
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FIG. 10: The factor Cog = x/^^ appearing in equations ([5| 
and ([8| in Section |Tl] Solid lines from left to right correspond 
to beams of energy 200 MeV, 1 GeV, and 6 GeV, respectively, 
incident on a Tungsten target. The dashed line corresponds 
to a 6 GeV beam incident on an Aluminum target. 



density (in g/crn^) and unit radiation length (in g/cm^) 
of the target material, respectively, and 



dP{t) 
di 



1 



(B2) 



is the differential decay probability where = 7ct 



Eict/itia' is the A! decay length given in (12 1. Also 

..bt-i bt 



I{Ei;Eo,t) 



SiE, - Eo 



T> 1 
T< 1 



(B3) 



is the energy distribution of electrons at position t in the 
target, where y = and b = 4/3. 

We can perform the t-integration explicitly in the limit 
of a very thin or thick target (T ^ 1 or T ^ 1). For a 
thin target, we find 



dNthin 

dxdz 



Nqp£q da 
A dx 



e - 1 



dP{z) 
dz 



(B4) 



For a thick target, we neglect the ^-dependence in the A' 
decay probability in (Bl I, since most production occurs 



within the first radiation length and thus well before the 
end of the dump. Here we find 



dNi 



thick 



dxdz 



Ne 



NoXo 



A 

Eq da 
Eydx' 



dEiI{Ei;Eo,T) 

dP{z) 

,_Eai dz 



(B5) 



where 



i{Ei;Eo,T) 



dtI{E,;Eo,t) 



1-l-y" (bTlny- 1) 
Eoby{lny)^ 



Eoby(liiyy 



(B6) 



16 



as T Qo. For finite T, the Hmiting form is a good ap- 
proximation for small and moderate y, i.e. for electrons 
that carry a large fraction of the initial beam energy (for 
y < 0.5 and T > 7 it is correct to within 1%). 

We note that by ^ and (Bll, -^-^^ is proportional 

to 



1 Sa^'e^Z^x NqXo 



to 



A 



(B7) 



However, -rS- 



^rad — f{Z)) 

the ato 

of the target atoms |41], and we obtain 



Lrad, L'^ad^ / ^ogs Set by the atomic form factors 



dN 
dxdz 



min(r, 1) 



(B8) 



with only logarithmic dependence on the target nucleus 
Z. This expression has a simple physical interpretation: 
an electron is slowed in a radiation length by radiating 
a small number of relatively hard but coUinear photons. 
The probability of instead radiating an A' is suppressed 
by the squared ratio of the couplings, e^, and the squared 

2 

ratio of the masses because it requires a higher 

invariant-mass intermediate state. 

For decays at a given z, we applied an acceptance based 
on the detector geometry and energy cuts before integrat- 
ing to obtain the expected total rate. 



APPENDIX C: KINEMATICS OF SIGNAL AND 
CONTROLLING BACKGROUNDS 

The dominant QED backgrounds for A' production are 
the trident reactions shown in Figure |3] The 7* process 
contributes an irreducible background to A' . The 

Bethe-Heitler process has a much larger rate, but can 
be controlled by exploiting its very different kinematics 
compared to the signal. Our aim in this appendix is to 
quantitatively describe the singularity structure of the 
Bethe-Heitler process and derive an effective set of cuts 
on lab-frame observables. Of course a more accurate ba- 
sis for a final design would rely on monte carlo for these 
processes, but the simple calculation clarifies the origin 
of the large Bethe-Heitler cross-section and how to regu- 
late it by cutting away from the dangerous "forward" and 
"asymmetric" regions of phase space in the lab frame. 

As an important reference, we will start by recall- 
ing the kinematic properties of A' production and de- 
cay using the results of Appendix |Xj We again consider 
a monochromatic incident electron beam of energy E^). 
Let 6cm be the emission angle of the forward decay prod- 
uct relative to the A' direction in the A' rest frame. 
Let 0^1 be the emission angle of the A' relative to the 
beam direction in the lab frame. As we have shown, the 
characteristic A' emission angle is small and is set by 



where x = Ea'/Eq. In the limit of 



small the lab frame opening angles d± and energies 
E± of the A' decay products are, 



xE 

E± = ^(l±COs(0em)), 



(CI) 



tan(g±) = ±^-J ]ll2?r\ +tM0A'), (C2) 

7 y 1 ±C0S(6'cmj 

where 7 — . 

The characteristic transverse momentum of the A' is 
Pa'.± ~ EaiOai ~ ^J\ — xTTiAi , while the typical recoil 

of the target is Igmml ~ ( 2xEo ) "^^^ median value 

ofl — xis 1 — a; « ^^^^ for mA' ^ 50MeV, implying 
that |gmm| is parametrically smaller than pa',± by ~ 
\p^^- Evidently, the recoiling electron largely balances 

the recoil of the A' . The energy E'^' and angle On of 
the recoiling final state beam electron in the laboratory 
frame is. 



Er 
tan(6'i^) 



(1 - x)Eo W TTIA' 

ruA' 



rriA' 

En 



(1 



2En 



+ ...). 



(C3) 
(C4) 



Note the relatively wide angle of the recoilin g el ec tron 
relative to the A' decay products. Equations (C1)-(C4| 



summarize the important kinematic characteristics of A 
production. 

As with A' production, trident reactions can also be 
analyzed using the Weizsacker- Williams approximation, 
where we group an outgoing I'^i" pair with fixed invari- 
ant mass to act as the A' candidate. In the case 
of Bethe-Heitler production, we can further approximate 
the beam electron as splitting at small angle into the 
recoil electron and a nearly on-shcU photon, which scat- 
ters with the Coulomb photon into the £^£^ pair (see 
Figure 111. There are two sources of large logs in the 
Bethe-Heitler cross-section: the soft and coUinear logs 
in the photon radiation by the electron and a forward- 
scattering log in the 77 £~^£~ subprocess. Both are 
regulated by We, and can be further suppressed by kine- 
matic requirements. We postpone the derivation — a 
matter of standard results and kinematic bookkeeping 
— to the following section, and focus here on its impli- 
cations. 

In this approximation, upon integrating over an invari- 
ant mass window of size Sm about niA' , we find 



da 



2a^X Sm 1 1 + (1 - x) = 



dxd cos BA'dc 



\ r X 



1 + c 1 



1 



1 + c 



(C5) 



where x and Qa' are defined as before, and c = cos 9cm, 
where 6cm is measured relative to the axis of the incoming 
photons in the 77 — > £~^£~ process, which is near enough 
to the beam axis for our purposes. This displays the 
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expected singularities at small x (soft), |c| ^ 1 (forward 
scattering), and small Oa' (coUinear). This is quite differ- 
ent behavior from the signal and radiative backgrounds, 
which are peaked at large values of x and slowly varying 
in c. 

Requiring x above 1 — 5, with 5 near or below its me- 
dian value 5 = max(TO^/ / Eq, m^/mA') keeps a large frac- 
tion of the signal and suppresses the Bethe-Heitler back- 
ground by a factor of 5. Likewise, the signal is relatively 
flat in c = cos 0cm- According to (CI) and (C2|, we 



can ensure modest c by constraining the ratio of the lab- 
frame energies or of the opening angles of the two decay 
products to be near unity. 

There remains the coUinear singularity at small Oa'- 
The signal is also peaked forward, but with a singularity 
regulated by m^' rather than Wg, so that it is produced 
at much wider characteristic A' angles than Bethe-Heitler 
processes. Since the absolute angles ^ {niA' I Eof"^"^ ^'"'^ 
still small and must be obtained by summing two mo- 
menta, it is impractical to place a lower bound on the 
A! emission angle. It is probably much easier to cut on 
the angle of the recoiling electron. Since the recoiling 
electron has much lower energy at large x and approxi- 
mately balances the transverse momentum of the A', it 
is emitted at a much larger angle. For the median x, the 
electron energy ~ rtiA' and its angle ^ (to^//£'o)^/^. Re- 
quiring the recoiling electron momentum and angle near 
these values significantly reduces the Bethe-Heitler rate. 
In a similar spirit, it may or may not also be easier to 
implement a tight x cut using the energy of the recoil 
electron rather than the total energy of the decay prod- 
ucts making up the A candidate. After these cuts, we 
find 



IGa^X Sm 1 



TT niA' tUai 

X log ((1 - a;mm)"^^fl,mi„) (1 - X„^^n) 



(tanh 



./2) 



(C6) 



When Ocut is small relative to the typical angular spread 
-g^(l — x) of the signal and radiative backgrounds, the 
factor on the first line is related to the lepton-pair cross- 
section fJ^ad with the same cuts from the radiative dia- 
grams alone by 



QcTrad log 



^radi^R^rnim \c\max), (C7) 



where 1 — x,nax = max(TO^ /i?Q, m^/m^) is the value of 
X where the log divergence in arad is regulated, and e^ad 
the efficiency for the radiative process (or signal) to pass 
On and |c| selections. 

Let us consider a representative case, Eq ~ 5 GeV, 
niA' = 0.5 GeV. Requiring x > 0.9, Or > 1/10, and 
E+/E^ < 3 (i.e. c < 0.5), we retain roughly 20% of the 
signal and reduce the contribution of Bethe-Heitler to 
the signal region to roughly the same size as the radia- 
tive contribution. A proper optimization of these cuts 
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FIG. 11: Bethe-Heitler reactions viewed as hard 77 ^ i+r 
processes. 



is best done with full monte carlo for the background, 
and of course depends on the characteristics of an in- 
dividual experiment, but we have confirmed numerically 
that the kinematic differences between Bethe-Heitler and 
radiative production are sufficient that it can be made 
sub-dominant while maintaining high efficiency for the 
A' signal. 

In the case of A' e+e^, the Bcthc-Hcitler process 
can also contribute with the electron labeled £~ in Figure 
[s] identified as the recoiler and ^+ and e~ forming the A' 
candidate. Here the recoil electron kinematics is as in 
the signal process but there is a forward-scattering sin- 
gularity when most of the A'-candidate energy is carried 
by the e~ with a softer e^. The cuts above remove this 
singularity as well, and we will not discuss it further. 



Bethe-Heitler Pair Production in the CoUinear 
Splitting Approximation 

We now extend the earlier Weizsacker- Williams treat- 
ment to compute the Bethe-Heitler pair production cross- 
section. Using the notation defined in Appendix but 
now with an outgoing lepton/anti-lepton pair with mo- 
menta /l^, the fully differential cross section is. 



dajpi +Pi^p2 + l+ + l- +Pf _ /ax\ 
dE A' d cos 9 A' dm'^did(f>d V tt / 

da{pi+q^P2 + l+ + 1^) 



EqxPa' 
(l-x) 



d{pi ■ a)dm'^dtd4>d 



(C8) 



t=t„ 



where a = l'^ + I 



is the total four-momentum of the 
- their invariant mass, and t = 
(Z+ - g)^ = (/+ - P/ + Pi)2. 0rf the angle between pi 
and in the rest frame of a. As before, the hard sub- 
processes (2 ^ 3 in this case) are separated from the soft 
Coulomb photon exchange. 

Starting from equation ( C8 1 , we will analyze the be- 



havior of Bethe-Heitler reactions (see Figure [3] (b)) rela- 
tive to A! production in a leading logarithm approxima- 
tion. This suffices to identify the singularities. 

To write the right-hand side of ( C8 1 in the approx- 



imation of near-coUinear splitting, it is useful to intro- 
duce the more familiar z and pj^ variables and relate 
them to the kinematic variables in the center-of-mass 
frame of the beam electron and Coulomb photon and 
to the lab-frame. We recall that in the lab frame, the 
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Coulomb photon is purely spacelike, with \q\ = y/tmin ~ 
2Eo(i-'x) ~ m'^/{2E{jx) for small 9a' (we have replaced 
m\, with the invariant mass rn? for the ofF-shell pro- 
cesses, but the kinematic conditions are unchanged). As 
\q\ at tmm kinematics depends on x and 9a', the CM 
frame does too — it is obtained from the CM frame by 
boosting with (3 ~ 1 — q/ Eq (we drop much smaller 9^'- 
dependent corrections) , leading to CM-frame momenta 
for the incoming electron and Coulomb photon. 



where p 
tum 



(P,0,0,p) 
(/3p,0,0,-p) 



(C9) 
(CIO) 



T^TT^^^vcn- The recoil electron has momen- 



where 



dcr(77 l^l ) 
dtd4>d 



a" I t 

T4 



(C15) 



We emphasize that the recoil electron is always right- 
moving, so the splitting approximation is valid where 
the Bethe-Heitler cross-section is largest, but is never 
a good approximation for the radiative (or signal) pro- 
cesses, where the 7* e~e~^* "splitting" always pro- 
duces an electron going backwards relative to the 7* . 



by 



The kinematic variables of (C8) are related to z,p± 
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P2 



il-x)EQ,Eo9A'X,0 
(1 - x)Eo 



EoO%x^\^ (CU) 



2(1 -x) 

(zp,Pi,o,y^zV-Pi), (C12) 



= = (l + /3-2z)(H-/3)p2 



Pi-a 



<t>d = (l>d, 



2-2 p2 
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(C17) 
(C18) 



where 



z — I — X, p± = E{)9a'X, 



(C13) 



leading to a Jacobian factor 



up to corrections of higher order in 9^' and m^. For small 
p_L, we can now approximate 



du 



dzdp\did4>d 



\2'k) (1 



[1 + z^] dcr(77^?+r) 
(1 - z)p']_ did<j)d 



d{pi.a)diin? p'^{\ + (3) dp'^dz 



2x(l ~ a;) d^ 



(C14) Taking I + (3 ^ 2, p^ 
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4x 



we find 



dp-^dz 



(C19) 



da{pi + q ^ P2 + 1+ + 1-) 



d(pi ■ a)dm'^did(f>d 



a 1 + (1 - x)2 2x(l -x) ft u 
2^ El9\,x^ ^ ^ yii^ 1 



and hence by ( C8 ) , 



da 



a^X l + {l-xY I ( i u 



dxd cos 9 A' dm? dtdipd t^"^ (^a'X \u t 



(C20) 



(C21) 



Changing variables from t to c = cos 9 cm of the 2 — > 2 m? from m\, — niA'Sm to m\i + niA'Sm, we obtain the 
process, i = ^(1 — c), and integrating over (f>d and over result of Equation (C5). 
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